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Abstract 

In this paper a new functional integral representation for classical dy- 
namics is introduced. It is achieved by rewriting the Liouville picture in 
terms of bosonic creation-annihilation operators and utilizing the stan- 
dard derivation of functional integrals for dynamical quantities in the 
coherent states representation. This results in a new class of functional 
integrals which are exactly solvable and can be found explicitly when the 
underlying classical systems are integrable. 
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1 Introduction 

Functional integrals are very popular in quantum mechanics and quantum statis- 
tical physics as a basic tool of both perturbative and non-perturbative analysis. 
The most fruitful applications of the functional integrals include Faddeev-Popov 
ghosts in the field theory non-perturbative corrections to the standard per- 
turbation theory series [HE], functional formalism for disordered systems 0], 
introduction of collective variables based on Stratonovich-Hubbard transfor- 
mation(see and refs therein) and applications to algebraic and differential 
geometry 013 EI. 

The beauty of the functional integrals lies in the fact that they essentially 
reduce quantum mechanical problems to a classical mechanical picture with an 
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additional stochastic noise. This means that instead of necessity of deahng with 
operators and partial differential equations (sometimes infinitely dimensional) 
one can think about the problem in terms of classical observables and, at the 
final stage, to average them across possible trajectories generated by (quantum) 
noise. This observation allows to find a formal expression for any quantum 
observable and is a cornerstone of all mentioned above applications. 

If functional integrals bring quantum mechanics back to classics, it is possible 
to formulate classical dynamics in terms of quantum one and therefore benefit 
from a quantum-mechanical toolbox. The question is not new. In last ten years 
series of papers examined one of possible ways to introduce the path inte- 
gral formalism into the classical mechanics. Although bringing about interesting 
mathematics they, in our view, did not result in a workable quantum-mechanical 
formalism which would allow a user to employ many useful functional-integral 
tricks which proved to be so helpful in the theory of many interacting degrees 
of freedom, disordered systems and renormalisation group. In this paper we are 
trying to fill the gap and to introduce a functional integral representation of 
classical dynamics which is similar in many respects to the standard functional 
integrals of many-body theory and allows perturbative and non-perturbative 
treatment of classical many-body systems and classical random dynamics. As 
a by-product we also obtain class of quantum mechanical systems which are 
related in our formalism to classical mechanics and therefore allow exact solu- 
tions. 

The paper is organized as follows. In the next section we briefly remind the 
Liouville formulation of classical dynamics, outline the main idea of the paper 
and introduce the creation-annihilation pairs for classical momentum and co- 
ordinates. Section 2 treats in great detail the case of one-dimensional classical 
harmonic oscillator. It is shown that there exists a two-dimensional quadratic 
quantum mechanical system which is equivalent to the classical oscillator in the 
sense that all quantum dynamical charatecteristics can be expressed through 
the corresponding classical trajectories. Section 3 builds on this and describes 
a general n-dimensional construction in the case of arbitrary classical potential. 
This results in the description of a new class of functional integrals and the cor- 
responding quantum systems which are exactly solvable and which dynamical 
quantities can be found explicitly if the underlying classical systems are inte- 
grable. Section 4 concludes the paper with a list of possible application of the 
developed formalism. 
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2 Classical Mechanics in the coherent represen- 
tation 



Let us consider a mechanical system with the phase space Al — {p, q}^, algebra 
of observables A, set of states S, and Hamiltonian H. Assume that the phase 
space Ai of the mechanical system has a volume form invariant under the phase 
flow with Hamiltonian H. Denoting this volume form by dx, we can write a 
probability measure d/i as dfi{x) — p(x,t)dx, where p{x,t) is a positive distri- 
bution (generalized function) on A4. For the pure state supported a,t xq € M 
we have p{x,t) = 6(x — xq) Dirac (5-function. This representation introduces 
Liouvilles picture, in which observables do not depend on time 

1 = 0, feA, 
at 

and states dfi{x) = p{x, t)dx satisfy Liouvilles equation, 
dp{x,t) 



dt 



~ —{H, p}, p{x, t)dx £ S, 



which is understood in the distributional sense, and Poisson brackets is given 
by: 

dp dq dq dp 
We can rewrite this equation in Schrodinger manner form: 

« Q-^ ^ ^P,<l)P\P,q]t)■ (1) 

where C is Liouville operator determined as 

.dH_^ .dH_d_ 
dp dq dq dp' 

Every solution of the initial value problem: p{p,q',t)\t=Q = p(p',q',0) for 
the Equation can be written in terms of the evolution operator U{t) = 
Q-itc^ which is a unitary operator that relates the state functions of the system 
p{p, q; T) with p{p' , q'; 0). Operator U{t) is an integral operator with the kernel 
U{p, q\p'q') and we can write: 

p{p, q)^ [ U{p, q\p', q')p{p', q') dp' dq' . (3) 



^To keep formulae simple in this section we do not stress that both p and q can be n- 
dimensional vectors. 
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The next step is to introduce creation and annihilation operators ai, 02, a[ 



ak: 



ai = 



02 



d 



P- 



V2\ dq 
with standard commutation relations 



[ai,a\] = 1, [02, 4] = 1' 

all other commutators are 0. These operators allow us to define the corre- 
sponding coherent states which are eigenstates of the operators a = (01,02), 
at = (a|,4): 



\z) = \z1z2) = 



E 



n=0 



ml 



|0)|0)e- 



_ xn 00 / — \m 



(^l^l) V- ("2^2) l|^^|2 



(^1 = {z,z,\ = (0| (0| 5: 5: 

■^^ — ' n! -"^ — ' ml 

Here |0) is defined by the equation o |0) =0 and unit decomposition is: 

dzdz , , , , 



(4) 
(5) 



27ri 



z iz 



It is not difficult to find the coherent states in the coordinate representation by 
multiplying this equation from the left with {x\ we now get (a;| a |0) = 0, or 



1 

71 



x + 



d_ 
dx 



) = 0. 



(x| 0) 



This differential equation is simple to solve, and we obtain 

ix 1/4 

In similar manner one can find the a;- representation of other eigenstates \n) by 
applying the operator (a^)" to the groimd state: 

1 



{x\ n) = yx 

p^'l^ (x 



-6 2 



Here Hn{x) = e^'/^ (a; — e ^ is the n-th Hermite polynomial. It is now 
straitforward to write down the coherent states in a;-representation: 

2 \ 1 — l_g 2 

^n!<!/7r\/2^ 

Tl=0 ^ 

-+\/2xz( 



{X\ Z) 



. 1^1 \ ^ Z 



n=0 



-e 2 2 



n=0 



1 

e 2 2 6 2 



^(6) 
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where we used well known formula for the generating function of Hcrmite poly- 
nomials: X^J^o ^"(^)fn" ~ exp {— + 2sa;}. The same analysis yields the 
following expression for the conjugated coherent states: 

(z| x) = ^ exp _ ^ „ 5 + y2xzl (7) 



The coherent states are the basis for derivation of the functional integral. As 
soon as the Liouville operator |2Jl is expresses in terms the creation-annihilation 
operators ai,2, o\ 2 in the so-called Normal form (i.e. when all creation operators 
stand to the left of the annihilation operators), £„orm(a^,a), one can find the 
evolution operator U{p, q\p' , q') as 

U{p,q\p',q') = {x\UiT,0)W) 

dzdzdz' dz' {x\ z) {z\ U{T, 0) \z') {z'\ x') , (8) 



where dzdz = ^^^^ and dz'dz' = and the kernel of evolution 

operator is given by the following functional integral [5]: 

{z{z2\U{T,Q)\z'^z'2) - j I?Vi2??i2?^22?'/'2 
exp \j2MT)MT) - \\MT)? - 



,k=l 



T 



+ (-V'i(i)^iW-V'2WV'2(i)-*'C„o™(?Ai,i^2,Vi,^2))di>. (9) 



Here V-ip^Dipi — J| I?V'2^''/'2 — Y\ '^^'2^^^ and initial conditions are: 

t t 



V'lW^zj, ^2(0)- 4 
V'i(T) = zr, '4>2{T) = z^ 



This is the functional integral representation of the classical dynamics we study 
in the paper. 



3 Harmonic oscillator 

Before moving further, we would like to consider a simple example which became 
rather customary to test new techniques, the harmonic oscillator. In this section 
we apply the derived above functional integral representation for the evolution 
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operator of classical harmonic oscillator with hamiltonian H{p, q) = ^ (jp + q^)- 
The kernel of the classical evolution operator or propagator gives us solutions 
of the classical dynamical equations: 

9H . dH 

or 

U{q,p\q\p')^{q,p\q',p')^5iq-q')S{p-p'), T>0 (11) 

with 

p' = p' cos LuT — q' sin ojT, q' = q' cos ujT + p' sin luT. 

and thus completely describes system behavior. We now obtain these expres- 
sions from the functional integral for the evolution operator. 

Using Equation Q one can derive the following expression for the Liouville 
operator: 

C{p,q)^^u^q-^p-j , (12) 

which can be easily expressed in the terms creation and annihilation operators 
Ui, 0,2, fl| , 02: 

If d\ I 1 



1 / d\ t w 9 



as 



£(01,4,02,4) = — (^(4 + a2)(ai - 4) - (a| + ai)(a2 - 4)^ 

= (^4«l - 4«2) = £Aform(ai, 4j 12, 4)- (13) 

This operator is already in the Normal form and can be directly used in the 
functional integral for (z|Z^(T, 0)| z'): 



exp ij^ii.mMT) - ^|?fe(r)p - i|^fe(0)p 

+ / (-?^i(t)^l(t)-^2W^2(t)+C^(?2W^lW-V^lW^2(i)))rfai4) 
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where V^piVipi — Y[ I??/'2^V'2 = '^^■ziTi''' with the mitial conditions: 

t t 

The functional integral is Gaussian and can be easily evaluated by finding the 
saddle point trajectory which are defined by the equations: 

Solutions of these equations are: 



and 



ipi{t) = z'2 cos Lot — z'l siiiLot, 
■02 (t) — z[ COS ujt + Z2 sin Lot^ 



V'i(i) ^2 c.osLo{t - T) - z'l svauj{t - T), 
-02 (i) = ^1 cosw(t — T) + sinw(t — T). 



Therefore the integral turns into: 



(Z1Z2I ?7(r, 0) \z[z2) = exp <j (ziz^ + Z2Z2) cosljT - ^-^^ 



-(zi4 + Z2z;)sina;T^^-^^ (15) 



Combining together equations ©JO,® and (|15|) one can get the following 
expression for the kernel of the evolution operator: 



dzidzi dz2dz2 dz'i dz^ dz'2 dz'- 



exp { - |Z2|2 -\z[\^^\z'2\^^^-^-^-^-^-^-^-^ 

^' ' ' ' ' '1' 2 2 2 2 2 2 2 2 



exp 



|\/2(ziP + Z2q + z'j^p' + z'2q') + {ziz' + Z2z') cosloT — {z'iz'2 + ^1^2) sinajr| 



Now it is possible to show that this is Dirac 5-function on the solutions of 
the classical mechanical equations of motion. Indeed, consider the following 
decomposition: 
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5{x - x') = {x\ I \x') = / (x| z) (z| z') (7| x') 



4:im.TT 



I 



dzdzdz'dz' 1 f UP x"^ z^ 



cxp < --^^ — + V2xz > X 



4i7ri7r 2 2 2 

M|2 -2 _ |,|2 



exp < -^^T — + V2x'z' 



2 2 2 2 2 J 

Using this decomposition and introducing a set of new variables: 

z[ = z[ cos ujT — ^2 sin luT, 'z^ = z'^ cos loT — zJ, sin loT, 
Z2 = Z2 COS (jjT + z[ sin wT, 'z-^ = z'l cos wT + Z2 sin wT, 

and 

^ = p' cos wT — sin ivT, = q' cos + p' sin wT. 
the kernel transforms into 

odZo 
^ ^ X 



, , ,s A dzidzidz2dz2dz{dz-,dz'o 
U{p,q\p,q) = J ^5 

exp{-|zi|2-|z2r-|£ip-|z~2'P} 



X 



^2 2 /2 /2 2 2 jy2 ;j/2 
exD < _£2_£2__£2__P__Q^_P ^>x 

^'22 2 2 22 2 2' 

exp I \/2(2;ip + ^;2g + "zi^ + %q') + {ziz[ + Z2z'2) | . 

which produces the required result 

U{p,q\p',q')=6iq-q')5{P-p') ■ 

This shows that the functional integral representation for the matrix elements of 

classical evolution operator indeed reproduces the Liouville flow on the classical 
phase space and is equivalent to the solution of classical Hamilton equations. 

4 New class of localization formulas for func- 
tional integrals 

In this section we describe a new class of functional integrals which can be 
evaluated exactly using their relationship to some classical dynamical systems. 
The main result is summarized in the following 
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Theorem 1 For any function V{qi . . . g„) the following integral: 

Iv{z,z)= [ P^iPV^iP^22?V'2exp|f^^f,(T)^UT^)-^l?fe(T^)l'-^IV'U0)|4 

•' U=ifc=i J 

' n 

Y^-i^,mi{t)-ipi{t)r2{t) + [i^,{tmt) -i^lmm) 

(Vi(t)-?i(i))4S™(V'>V^)}(16) 



exp ■ 



T n 

+ 



where VtpiVipi = I?V2^^2 — Yl '^^27!^^ ^^^'^ initial conditions: 



and Ci^norm{a, at) = AT ^ |X - 
is equal to 



9' = 75(4+4') 



g2 2 22g2 2 2 2 



'i2 — i2 — i2 



e 2 i 1 |-gV:^^^l P -+^2 9 U )+^2r(i )) j ^ (^ij^ 

where {p*(T), g*(T)} are the solution of the classical Hamiltonian equations with 
the Hamilton function 



1 " 

H=-Y,P'' + V{q'...q'') . 



for the time t = T with initial conditions at t = being {p'*, g'*}. 

Proof. Let us consider classical system with the following Hamilton function: 

i—l i—1 i=l 

The corresponding Liouvillc operator C will have the form: 

Expressing this operator in terms of creation- annihilation operators a\,a2, , 02 



9 



one can get 

dV 



j=l * i=l ^ 



f dV_ 



9' = 7f(4+4') 



1 yoqi 

1=1 

If the this equation is rewritten in the Normal form: 

n . n 

j=i ^ 1=1 

n 

V2 



(21) 



(22) 



where 



9' = :^(a^+a+') 



it is straightforward to put down the functional integral for the matrix elements 
of the evolution operator in the coherent state representation: 

{ziZ2\U{T,Q)\z[z'^) = j PVi2??i2??22?V'2 
U=l fe=l J 



- {ij;\{t) - ?i(t))A,„(V', V^)'ii} (23) 

n —I I n —I I 

where ©Vi^^V-i = 11 11 ^^^^^ 2?V'22^V'2 = 11 11 with the initial condi- 

tions: 

M^) = A. ^2(0) = 
V'i(T) = zr, '02(1')=^ 

The same matrix element can be found as 

U{{A,A]\{z'l,z'^]) = {{z{,zi} \U{0,T)\ {z'^,z','}) 

Yl dp'di/'d<fdq'^ {{z\,zl} I {p\ q'}) {{p'\ q"} \ {z^ 4}) X 
i=l 

({p\g*}|W(0,T)|{p^gn) • (24) 



10 



Using equations © and iQ one can easily get: 

lldp^dp'^dq^dq'^ {{p\q^} \U{0,T)\ {p'\q''}) X 

TZ {{p\ q^}, {p'\ q'^}, {4, 4), {zi, z^}) , (25) 



where 



7^ g^, 9"}, 4}, ^2}) 
^{{z\,z\}\{p\q^}){{p'\q'^}\{4,4}) 



n 



2 ivi|2 



1 _p^_g^_p^_^ 

g 222 2g2 2 2 



^/i2 2/i2 — i2 — i2 

g T T J i-gV2(^rp"+22'9"+^lP'+^2'3' 



(26) 



Taking into account that the classical evolution operator is a delta-function on 
the solution of the classical equations it is easy to see that expression (|25|l can 
be re-arranged as: 

g-^-^-^-4^p-4^-4^-¥-¥pV2(z;>"+4>"+-lp'(T)+-29'(T)) 



where {p*(r), 9'(T)} are the solution of the classical Hamiltonian equations for 
the time t — T with initial conditions at t = being {p'*,g"}. This gives 
the formula for exact evaluation of the described class of functional integrals 
in terms of solutions of the corresponding classical equations of motion. In the 
case of integrable systems when the solutions can be written down explicitly, 
the functional integrals will be evaluated explicitly as well. This completes the 
proof. 



5 Concluding Remarks 

In this paper we suggested a new functional integral representation for classi- 
cal evolution operator. Having calculated this expression one can study time 
evolution of any classical observable as well as statistical properties of classical 
systems. As a by-product we also found a set of quantum mechanical systems 
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which are related to the classical mechanics and can be solved exactly by reduc- 
ing the problem to solutions of the corresponding classical equations. 

As we pointed in the Introduction the existence of the functional integral 
representation opens interesting technical opportunities. Let us list here the 
main directions which, we believe, promise interesting results in future: 

1. Study of classical dynamical systems in the presence of random forces. 
Since the evolution operator can be formally written for any interaction 
structure, one can derive directly a functional integral representation for 
averaged dynamics of classical observables as well as their statistical char- 
acteristics in close analogy with the corresponding quantum case 4 . One 
of applications in this case would be study of classical chaotic systems by 
means of non-perturbative field-theoretical methods. 

2. Collective variables in classical dynamical systems. One can apply func- 
tional form of Stratonovich-Hubbard transformation to introduce collec- 
tive variables into the classical systems with many interacting degrees of 
freedom. This might be an interesting bridge between classical n-particle 
dynamics and classical statistical physics. 

3. Investigation of possibility to introduce dissipative forces into the for- 
malism. This can potentially lead to a new instrument of study non- 
equilibrium open systems and self-organization. 

4. Exact solutions of quantum mechanical or field theoretical problems in 
terms of the related classical trajectories. One can see the relation estab- 
lished here between as certain class of quantum systems and some classical 
systems as a form of localization formulas for the corresponding quantum 
functional integrals. It is interesting to see how wide is the class and if 
the formalism can be generalized to include even more quantum systems 
(not necessary 2n-dimensional ones). 

We hope to return to these issues in our future publications. 
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